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Abstract 

We consider the polymer quantization of the Einstein wormhole throat theory 
for an eternal Schwarzschild black hole. We numerically solve the difference 
equation describing the quantum evolution of an initially Gaussian, semi-classical 
wave packet. As expected from previous work on loop quantum cosmology, the 
wave packet remains semi-classical until it nears the classical singularity at which 
point it enters a quantum regime in which the fluctuations become large. The 
expectation value of the radius reaches a minimum as the wave packet is reflected 
from the origin and emerges to form a near Gaussian but asymmetrical semi- 
classical state at late times. The value of the minimum depends in a non-trivial 
way on the initial mass/energy of the pulse, its width and the polymerization 
scale. For wave packets that are sufficiently narrow near the bounce, the semi- 
classical bounce radius is obtained. Although the numerics become difficult to 
control in this limit, we argue that for pulses of finite width the bounce persists 
as the polymerization scale goes to zero, suggesting that in this model the loop 
quantum gravity effects mimicked by polymer quantization do not play a crucial 
role in the quantum bounce. 
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1 Introduction 



There exists a large body of evidence supporting the notion that loop quantum 
gravity has the potential to resolve the singularities that plague classical cosmol- 
ogy ED El El 111 E] and black hole physics EJ El El EH EE] • Most of this work uses 
the polymer representation of quantum mechanics [121 EI] as a way of mimicking 
the features of a more complete theory of loop quantum gravity. It is important to 
note that polymer quantization is a viable quantization scheme in its own right, 
one that is unitarily inequivalent to the usual Schrodinger quantization. More- 
over, one of its defining features is a fundamental discreteness at some intrinsic 
microscopic length scale. It is plausible that an underlying discreteness at short 
distance scales is a general feature of any theory of quantum gravity, including 
loop quantum gravity, so that the qualitative features that emerge from poly- 
mer quantization may give clues regarding the short distance effects of quantum 
gravity, irrespective of its detailed microscopic structure. In this spirit, poly- 
mer quantization has been applied to a variety of different quantum mechanical 
systems including the harmonic oscillator [12] . Coulomb potential [H] and 1/r 2 
potential [15]. An important question that arises in the context of singular po- 
tentials concerns the extent to which polymer quantization provides a mechanism 
for singularity resolution that is fundamentally distinct from Schrodinger quan- 
tization. 

We focus on a model for the interior of an eternal Schwarzschild black hole 
that was first developed in the quantum setting by Louko and Makela [16]. The 
dynamics is that of the minimum radius of the throat of the Einstein- Rosen bridge 
as a function of the proper time of a comoving observer at the throat. This model 
has two advantages over other models of black hole interiors that have appeared 
in the recent literature. First, it is fully reduced and describes the evolution of 
an eternal Schwarzschild black hole in terms of geometrical invariants. Second, 
it produces relatively simple equations at the classical and quantum level that 
are amenable to both analytic and numerical analysis. One disadvantage of the 
model is that it is difficult to reconstruct the full black hole quantum corrected 
spacetime from the fully reduced model, as was done in ITT] . 
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The current work is a continuation of [T7] in which the full polymer quantized 
throat theory was constructed and the energy spectrum was obtained using a nu- 
merical shooting method. In the present paper we numerically evolve an initially 
Gaussian wave packet in the polymer quantized throat theory. The initial packet 
is chosen to be semi-classical in the sense that its width is large compared to the 
polymer scale but small compared to the only other scale in the problem, namely 
the ADM mass. To the best of our knowledge, the analogous calculation for this 
model using Schrodinger quantization has not been done but our results can be 
compared to the semi-classical limit of the Schrodinger throat theory obtained in 

Ref. ma. 

The paper is organized as follows: In Section 2 we review the throat theory 
that we use to describe the classical dynamics of the black hole interior. Section 
3 summarizes previous results on the effective polymer dynamics, which yield the 
semi-classical spectrum and time evolution of the throat. Section 4 describes the 
full polymer theory as well as our numerical methods and results. Section 5 closes 
with a summary, conclusions and prospects for future work. 



2 Classical theory 



We consider a Hamiltonian system with a two-dimensional phase space and the 
Hamiltonian 

H =l(- + r), (2.1) 



2 V r 

where the configuration variable r takes positive values and p is the conjugate 
momentum. The equations of motion reduce to 

r> = ^ - 1, (2.2) 
r 

where M is the conserved value of H and the overdot denotes derivative with 
respect to the time t. Note that M is necessarily positive. 



The equation of motion (2.2) is identical to the equation of a radial timelike 
geodesic that passes through the bifurcation two-sphere on a Kruskal manifold of 
mass M, provided r is identified with the area-radius of the two-sphere and t is 
identified with the proper time on the geodesic [TS1 H2]- This suggests that the 
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Hamiltonian (2.1) could be interpreted as the Hamiltonian of the Einstein- Rosen 



wormhole throat. It was shown in [16] that this interpretation follows from a 



direct derivation of the Hamiltonian (2.1) from the geometrodynamical action of 



general relativity, after reducing the theory to spherical symmetry and anchor- 
ing the evolution of the spacelike hypersurfaces into the two asymptotically flat 
spacelike infinities of the Kruskal manifold in a particular fashion [20] . To ensure 
that the value of the Hamiltonian equals the Schwarzschild mass, the evolution at 
one infinity is taken to proceed at unit rate in the asymptotic Minkowski time t, 
while the evolution at the other end is frozen. To express the Hamiltonian in 
terms of the wormhole throat radius, the hypersurfaces are chosen to coincide 
with those of the Einstein- Rosen wormhole near the wormhole throat, and the 
time evolution is chosen so that t coincides with the proper time at the wormhole 



throat. Under these conditions the Hamiltonian duly takes the form (2.1), with 
the variable r being the area-radius of the wormhole throat [16] . A sketch of the 
foliation is shown in Figure [T] 

To polymer quantize the theory, a choice must be made for the variable to be 
polymerized, and this choice can affect the qualitative properties of the quantum 
theory [21 [3j. For reasons described in [T7] we work with the canonical chart 
(4>, n), where 



(2.3a) 

n = \l , (2.3b) 

with < cf) < oo and — oo < n < oo. The Hamiltonian takes the form 

H=^- (4n 2 + 1) . (2.4) 

We work in Planck units, G = c = h = 1. The conversion to geometric units, in 
which G — c — 1 but [h] = (length) 2 , is discussed in [T7] . 

3 Effective polymer theory 

As discussed in more detail in ET] . a semiclassical, or effective, polymer ap- 
proximation is obtained by replacing n i— > sin(/xll) ///in the classical Hamiltonian, 



3 



with the result 



H, 



a/2 



eff 



4sin 2 (/TI) 



+ 1 



(3-1) 



where fx is the polymerization scale. The resulting equations of motion are 
thought to approximate the time evolution of the expectation value of x and 
p for semi-classical states. These equations read 

10 1 / 2 sin(fiU) cos(>n) 



= {0, #cff} 



tl = {U,H cS } = - 



M 



where M denotes the constant value of H efi on a solution. 
It follows that has two turning points where = 0: 



2M 
2M 



Solving for in terms of 0, fx and M gives 



2M 
0V2 " 

(4 + fx 2 ) 



4 

,1/2 



+ 



2M 

0V2 ~ 



■1/2,1/2 



By integration we obtain 
2 



(2M-0V2)(0i/2 -2Mk) 



/2M(1 + k)- 
+ M(l + k) arcsin I - 



a/2 



+ C, 



(3.2) 
(3.3) 



(3.4) 
(3.5) 



(3.6) 



(3.7) 



2M(1 - k) 

where k = fx 2 / (4 + fx 2 ). This equation gives implicitly the trajectory 0(t). Equa- 
tions (3.4), ( |3.5[ ) and (3.7) are the key ingredients from the effective theory that 
we will need for comparison with the time evolution in the full polymer theory. 



4 Full polymer theory 

Here we summarize only the features of the quantum theory that are needed for 
the subsequent calculation. We refer the reader to [121 [131 EH] and references 
therein for further details. 
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4.1 Operators 



We take the polymer Hilbert space % to have the orthonormal basis {\mfi) m G Z} 
of eigenstates of the operator 0, 



4>\mfJ,) = rnji \rnfj,), 



where /i is the polymerization scale as in Section [3] The orthonormality relation 
reads 

(mfi\m'fi) = 5 mjm t. (4.2) 

The momentum operator IT is defined by 

1 



where tin is the finite translation operator, 



Un\mn) = \ (m+ l)n). 



(4.3) 



(4.4) 



We note in passing that equations (4.3) and (4.4) correct a sign error in the 



definition of II in [T7j (see [T2] ) ; this error does however not affect the conclusions 



in [T7] as the polymer Hamiltonian (equation (4.5) below) is quadratic in II. 
We order the polymer Hamiltonian operator as [T7] 



^ pol = i(4fl0 1 / 2 fl + 1 / 2 ). 



(4.5) 



As is not positive definite, the square roots in (4.5) are not a priori defined. 
Following [T7], we define these square roots as 



1//2 |m/x) := |m/i| 1//2 |m/i). 



(4.6) 



The action of H po \ on the basis states then reads 



H vo i\mjj) 



2/i 3 / 2 



[\m + 1| 1/2 + |m - 1| 1/2 + /i 2 |m| 1/2 ) \mfi) 
-\m+ l| 1/2 |(m + 2)/i) - \m - l| 1/2 |(m - 2)/x) 



To solve Schrodinger's equation, 



d 



i-V(t) = H po Mt), 



(4.7) 



(4.8) 



we write 

- oo 

*( £ )=r 2 [/m(*) +*</«,(*)] l™/*>, (4-9) 



m=— oo 



c 

where / m and g m are real-valued functions of t and C is a real-valued normal- 



ization constant. Eq. (4.8) then reduces into a system of coupled first-order 



differential equations for the coefficient functions {f m } and {g m }. 

4.2 Initial state 

We take the initial state to be 

1 oo 

#(0) = |0 O)7 r o ) = 1 e -5 2 (m^o) 2 /(2M 2 ) e i TO | m/i ^ ( 410 ) 



c 

m=— oo 



where C is the normalization constant and the parameter 8 controls the width of 
the state. Using the notation above, 

j oo 

*(0)=~ lfm{0)+ig m (0)]\mfi), (4.11) 



C 

m=—oo 



with the identifications 



/ m (0) = e -^-A>)W) sin(7r m/i), (4.12) 
9m (0) = e -s 2 ^-M 2 /(2^) cos ( vromyU ) ; (4.13) 

This state is a generalization of a coherent state of quantized harmonic oscillator: 
it is peaked at the classical phase space point (0 o ,7r o ) p21 [21]. Note that in the 
present context the state can only be considered semi-classical when its width is 
large compared to the lattice spacing \i. To leading order, 

(A*) = ^ (4.14) 

where A0 is the standard deviation of 0. Hence the condition for the state being 
semiclassical is that y/2S -C 1. 

We start the time evolution when (0) is near the outer turning point, i.e. the 
horizon, which requires that 7r ~ 0. Numerically this means making tt so small 
that its further decrease does not change the results up to a desired accuracy. 
The relevant initial parameters are therefore </>o and 5. 
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4.3 Expectation values 



We have set up the polymer quantum theory so that the spectrum of <j) consists 
of all integer multiples of the polymer scale fi, including the negative integer 



multiples. This enabled us to define the momentum operator II (4.3) and the 



Hamiltonian operator H po \ (4.5) in terms of the translation operator (4.4). As 
H po i commutes with the parity operator, |m//) i-> \(—m)fj), the quantum theory 
decomposes into superselection sectors that consist respectively of symmetric and 
antisymmetric states under parity. As the classical theory satisfies < <p < oo, 
each of the two superselection sectors can be regarded as a distinct quantization 
of the classical theory: in the limit of small /i, numerical evidence indeed indicates 
that the spectra in the two sectors reduce respectively to those of the Schrodinger 
quantized theory with a corresponding boundary condition [T7] . The emergence of 
the symmetric and antisymmetric superselection sectors is qualitatively similar to 
what was found the Coulomb potential in [T3] and for the inverse square potential 

For computations in each of the superselection sectors, we may define the 
'physical' states ip(t) by excluding the negative values of m, 



I - 1 

ip(t) := (t) - - [/"»(*) + Wm(t)] \mn), 



(4.15) 



and we may normalize the inner product of the states tp^ and ip^ to take the 
form 



<^ (1) |^ (2) > 



1 



(/o (1) -^ J )(^ J +#) 



(i)\M 2 ) 



,(2)> 



m=l 

A normalized state then satisfies 

1 oo 



(4.16) 



(4.17) 



m=l 



1 Construction of polymer equivalents of more general Schrodinger boundary conditions is 
currently under investigation [55] . 
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The expectation value ((f)) and the variance A0 take the form 

(4>) = 



Q2 ^ > ^ ^ m ) ' 



m=l 



C 2 



.. oo 

^(0) 2 (/o 2 +^ 2 ) + E(^-^) 2 (^+^) 



m=l 



(4.18) 
(4.19) 



4.4 Numerical methods 



To obtain the coefficients f m and g m at a given time instant t, we solve the 



Schrodinger equation (4.8) using the fourth order Runge-Kutta method. Using 



Eq. (4.7) we obtain 



dfm , .dg 



(4.20) 



where we have arranged the coefficients f m and g m into the vectors / and g and 
we have defined the functional F m [/] as 



Fm[f\ 



2/x 3 / 2 _ 



^| m + l|l/2 + | m _ 1 |l/2 +/i 2| m |l/2j /) 

-Im+ll^/m+s-lm-ll 1 / 2 ^ 



(4.21) 



In the fourth order Runge-Kutta method, the time evolution of the coefficients 
f m and g m is given by the equations 

At 



f m (t + At) = f m (t) + =■ M m) + 2*T + 2 ^ mJ + k 



,M . , (m) 



where 



# m (t + At) = g m (t) 



Am) 
(m) 



6 L 
At 



(m) ,(m) 



S m ) 



F m [g] , 

Fm[g + "TplJ' 
„ r „ At ri 

F m [g + Atk 3 ], 



(4.22) 
(4.23) 



(4.24a) 
(4.24b) 

(4.24c) 
(4.24d) 
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and 



(4.25a) 



4 m) = ^/ + f^L 

ht ] = F m [f + ^h,l 
ht ] = F m [f + Ath 3 ]. 



(4.25b) 



(4.25d) 



(4.25c) 



To calculate the coefficients f m and g m at the next time step t + At, one therefore 
needs information from the neighboring lattice points at the time instant t. For 
instance, to obtain f m (t + At) one needs to know the elements of g(t) up to the 
lattice points m ± 8. This creates a boundary problem which is dealt with as 
follows. 

The boundary conditions at the origin pose no problem, as explained in the 
previous section. We can impose either symmetric or anti-symmetric boundary 
conditions and evaluate the relevant terms in the Runge-Kutta method for both 
positive and negative m. The results exhibited below used symmetric bound- 
ary conditions, but nothing substantial changes in the time evolution for anti- 
symmetric boundary conditions. The boundary at the large m limit is trickier. 
At every time step, we must extrapolate extra lattice points outside the original 
lattice in order to evolve f m and g m at the boundary. Extrapolating 8 new lat- 
tice points at every time step makes the code unstable so a different approach is 
needed. An easy and working solution is to use linear approximation 



for the lattice points m max — 7, . . . , m max . In doing so, one only needs to ex- 
trapolate two extra lattice points at the boundary which significantly reduces the 
numerical error. A simple linear extrapolation seems to be sufficient to obtain 
values for f m and g m on the two extra lattice points. 

The reliability of the code was confirmed by checking energy conservation and 
unitarity of the time evolution, as well as by the apparent agreement with the 



fm(t + At) 



f m (t) + AtF m [g] 



(4.26) 
(4.27) 
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semiclassical theory. To obtain a notion of energy describing the black hole mass, 
we calculate the expectation value of the Hamiltonian on the positive half line, 

1 



<M> 



C 2 2/z 3 / 2 



(/ 2 +^)-(l + S )(/0/2+^ 2 ) 

oo 

+ Yl Cfm + 9m) (y/m -1 + Vm + 1+ fi 2 ^n) 

9m9m+2 



(4.28) 

where the symmetry factor s = 1 for symmetric boundary conditions and s = —1 
for antisymmetric boundary conditions. The deviation of (M) from the classical 
value, M = a/0o/2, gives a quantitative measure of the extent to which the 
initial state is semiclassical. As one might expect, the initial Gaussian wave 



packet (4.10) is semiclassical as long as the width of the state is large compared 
to the polymer scale and the wave amplitude is close to zero at the origin. In our 
investigations, the deviation from the classical mass was typically of the order of 
0.5 %, although deviations up to 1.5 % were considered. In all these cases we 
found strong qualitative agreement with the effective polymer dynamics. 

Numerical accuracy was monitored during each run in three ways. First the 
expectation value of the Hamiltonian was calculated at each time step. For all 
runs the relative change in (M) was kept between 0.1 % and 0.15 %. Typically 
this required time steps of the order of 1 x 10~ 7 units of Planck time. Secondly, 
the norm of the state was calculated at each time step. Using the time steps sizes 
as mentioned above, the relative error was also within 0.1 % - 0.15 %, indicating 
the unitarity of the time evolution to the required order. Finally we checked on 
the convergence of the bounce value of (0) (denoted henceforth by (p m m) with 
decreased time step. The convergence of min in the limit t — > was found to 
be almost linear, which made the error estimate rather easy. For the given sizes 
of the time steps, the estimate for relative error in min was consistently below 
0.15 %, which is more than adequate for the purposes of our study. Smaller time 
steps would require significantly longer computation times and offer essentially 
no new information about the dynamics. 
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4.5 Results 



Figs. [2](5] show the time evolution of a Gaussian pulse with symmetric boundary 
conditions. The initial expectation value of the throat area for this run is 500 
/pj with an initial pulse width of 14 / P1 . The pulse initially spreads and then 
as it approaches the origin, it enters a quantum regime. As seen in Fig. [6] the 
expectation value reaches its minimum value of about 41 Planck units at time t 
38. The time evolution of (</>) is qualitatively in agreement with the semiclassical 



estimate obtained from Eq. (3.7). A notable detail of the quantum evolution is 
that the maximum at t « 76 is slightly below <f) = 500. 

The behavior of the width of the wave packet during the time evolution is 
somewhat unexpected. As seen from Fig. [7j narrow initial width of the pulse 
does not necessarily imply narrow width at the bounce. In particular there exists 
a minimum value for A0 at the bounce so that it cannot be made arbitrarily 
narrow by choosing the initial conditions. This appears to be a generic property 
for all initially Gaussian, semiclassical states. The width of the wave packet at 
the bounce (denoted henceforth by A0&) has significant effects on the nature of 
the bounce as described below. 

An interesting feature of the bounce is shown in Fig. |8j As the system ap- 
proaches quantum regime, the relative spread of the wave packet A0/ (0) increases 
strongly, and at the bounce region it becomes of the order of unity. For all the 
cases we were able to test, the maximum of A<p/ (0) was roughly between 0.4 and 
0.85, with the lowest values corresponding to the minimum of (cf. Fig. [7]). 
This means in particular that the uncertainty A0 becomes comparable to the 
expectation value of <fi at the bounce. Such strong fluctuations near the bounce 
naturally explain the relatively large values of the bounce radius that arise from 
the polymerized throat theory. It is also worth noticing that the relative spread 
of 4> makes a tiny dip shortly after the bounce. This feature seems to persist re- 
gardless of the initial values of the parameters 5, <fto and [i. Ultimately the reason 
for this asymmetric form of the curve stems from the finite width of the wave 
packet: the natural tendency to spread is balanced by the squeezing that occurs 
as the pulse moves towards the origin and occupies a smaller volume of configu- 
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ration space. This feature is not specific to the polymer theory but also arises in 
the Schrodinger quantization on the half-line. We have verified analytically that 
both free particle Gaussian wave packets and harmonic oscillator coherent states 
exhibit similar behavior for A0/0 as the packet scatters off the origin. 

The quantum bounce area m i n is generically somewhat higher than the semi- 
classical prediction. This also appears to be a consequence of the finite width of 
the wave packet at the bounce. Fig. [9] shows how </> min depends on the width 
of the pulse. As explained earlier, there is a limit of how narrow the pulse can 
be made at the bounce, and this, in turn, narrows the range of available data 
points. Nevertheless, it is quite clearly seen that in the limit where the width 
A(f> b gets small, min tends towards the value predicted by the effective theory, 



which for the initial parameters used in the figure equals 40 Z|j (cf. Eq. (3.5)). 
Since the relationship between min and Afa is almost linear, it is no surprise 
that Fig. [10] of min vs. A0 O is qualitatively similar to that of Fig. [7j Our results 
provide numerical confirmation of the range of validity of the effective theory 
in the semi-classical limit. For finite width packets, one obtains a minimum 
area that is determined by the width and is more or less independent of the the 
polymerization scale. This confirms the expectation, based on the study of the 
Schrodinger quantization of the throat theory [16J, that the central singularity in 
this model is also resolved in the continuum theory. 



The singularity avoidance in the continuum limit may be seen in Fig. 11 



which shows the minimum area m i n as a function of /x. The graph is obtained by 
keeping the physical width of the initial wave packet unchanged while decreasing 
the size of the lattice spacing /x. As /x becomes smaller, the number of lattice 
points increases, making the numerics increasingly more difficult to handle. Using 
0o = 80 and A0 O ~ 17.7 we were able to get down to /x = 0.075 with reasonable 
computation times. With the given initial conditions, we see that the magnitude 
of 0min does not change significantly when /x < 1. We were also able to confirm 
that the time evolution of (0) is qualitatively similar to that of Fig. [6] regardless 
of the size of /x. All this suggests that the existence of a finite polymer scale does 
not significantly alter the qualitative aspects of the singularity resolution in this 
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case. 



5 Conclusions 

We have considered the quantum dynamics of the Einstein- Rosen wormhole throat 
of an eternal Schwarzschild black hole using polymer quantization. An initially 
Gaussian, semi-classical wave packet remains semi-classical until it comes close 
to the singularity. The area expectation value reaches a positive minimum value 
as the packet is reflected from the singularity, and the packet emerges from the 
bounce to form a new, near-Gaussian semi-classical state. The bounce area de- 
pends in a non-trivial way on the initial energy and width of the pulse, and also 
on the polymerization scale: however, for wave packets that remain sufficiently 
narrow near the bounce, we find that the bounce area agrees with that of the 
semiclassical polymer theory |17j . 

Although the numerics become difficult to control as the polymerization scale 
decreases, our results provide evidence that for pulses of finite width the bounce 
persists as the polymerization scale keeps decreasing. We have not investigated 
numerically the behaviour of wave packets in the corresponding Schrodinger the- 
ory [TBI E]; but the unitarity of the Schrodinger theory implies that a bounce 
of some sort must occur, and we are not aware of reasons to expect that the 
small polymerization scale limit of the polymer bounce would differ from the 
Schrodinger bounce. 

For our wave packets, the wormhole area expectation value (0) is almost time- 
symmetric about the bounce, but the relative uncertainty A</>/((/>) is not, as seen 
in Figure |8j This asymmetry is a consequence of the spreading of the wave packet 
as it evolves. We have verified analytically that a similar asymmetry occurs also 
in Schrodinger quantization of a free point particle on the positive half-line with 
the corresponding boundary condition. In fact, for the free point particle the 
spreading is so strong that the packet must be carefully tuned to be narrow near 
the bounce in order to see a significant peaking in A0/ ((f)) about the bounce at 
all. 

Overall, our results do not suggest significant qualitative differences between 
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the polymer and Schrodinger singularity avoidance mechanisms for the Einstein- 
Rosen wormhole throat. In this respect the singularity of the throat theory 
appears to be no more insidious than that of the Coulomb potential [H] . We note 
that while the classical throat theory is completely deparametrized and involves 
no constraints, differences in polymer and continuum singularity resolution do 
arise in symmetry-reduced systems where constraints are present [H |2j [3j HJ EE EJ 
El 13 HOl [11] ; see in particular [4J for a situation where continuum quantization 
is argued to lead to physically unsatisfactory results. This suggests that the 
differences between polymer and Schrodinger quantizations may relate less to the 
singularities of the underlying classical theory than to whether the constraints 
are addressed prior to or after the polymer structure is introduced. 
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Figure 1: Schematic of the geometrodynamics of the time evolution of the throat 
of the Einstein-Rosen bridge as a function of proper time of a comoving observer 




Figure 2: Time evolution of wave packet, t = 0. In this figure fi = I, 5 = 0.05, 
0o = 500 and tt = -1 x 10" 7 « 0. 
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Figure 3: Time evolution of wave packet, t = 20. In this figure fi = 1, 5 = 0.05, 
0o = 500 and tt = -1 x 10~ 7 « 0. 
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Figure 4: Time evolution of wave packet, t = 35. In this figure /i = 1, 5 = 0.05, 
0o = 500 and tt = -1 x 10~ 7 w 0. 
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Figure 5: Time evolution of wave packet, t 
0o = 500 and vr = -1 x 10~ 7 « 0. 
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Figure 6: Evolution of ((f)) as a function of proper time t. The top curve (blue, 
solid) shows is for the full polymer theory whereas the bottom curve (red, dashed) 
is for the effective theory. The former generally produces a bounce radius that 
is slightly larger than the latter. In this figure fj, — 1, S — 0.05, 4>o = 500 and 
7T = -1 x 10- 7 w 0. 
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Figure 7: The value of the standard deviation of at the bounce, A0b, as a 
function of the initial value of the standard deviation A0o. in this figure \i = 1, 
0o = 1000 and vr = -1 x 10~ 7 » 0. 
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Figure 8: The value of A0/(0) as a function of the time. In this figure \i — 1, 
5 = 0.05, 0o = 500 and tt = -1 x 10" 7 « 0. 
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Figure 9: The value of the minimum of ((f)) as a function of For all the 

data points /i = 1, 0o = 1000, and n = — 1 x 1CT 7 ~ 0. Although cannot 
be made arbitrary small in the full quantum theory, the figure strongly suggests 
that formally in the limit A/fib — > 0, <f) m m approaches the value predicted by the 



effective theory, which for given values of parameters is 40 Z|j (cf. Eq. (3.5)). 




Figure 10: The value of the minimum of m ; n as a function of A0 O - For all the 
data points \i — 1, (f) Q — 1000, and 7r = — 1 x 10~ 7 w 0. 
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Figure 11: The value of m i n as a function of fi. All the data points correspond 
to the same the initial wave packet with a physical width A0 O ~ 17.7, <po = 80, 
and 7T = -1 x 1(T 7 w 0. 
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